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Abstract. In this paper we deal with the classical question of 
existence of polynomial in momenta integrals for geodesic flows 
on the 2-torus. For the quasi-linear system on coefficients of the 
polynomial integral we consider the region (so called elliptic re- 
gions) where there are complex-conjugate eigenvalues. We show 
that for quartic integrals in the other two eigenvalues are real and 
genuinely nonlinear. This observation together with the property 
of the system to be Rich (Semi-Hamiltonian) enables us to classify 
elliptic regions completely. The case of complex-conjugate eigen- 
values for the system corresponding to the integral of degree 3 is 
done similarly. These results show that if new integrable examples 
exist they could be found only within the region of Hyperbolicity 
of the quasi-linear system. 



1. Introduction 

Let p be a Riemannian metric on the 2-torus = M^/F, denotes 
the geodesic flow. Let : T*T^ be a function on the cotangent bundle 
which is homogeneous polynomial of degree n with respect to the fibre 
(notice that this condition is invariant with respect to the change of 
coordinates on the configuration space T^). We are looking for such 
an Fn which is an integral of motion for the geodesic flow p*, i.e. Fn o 

= Fn- This question leads immediately to a system of quasi-linear 
equations on coefficients of Fn and this is the aim of the present paper 
to study it for the degrees n = 3, 4. Let us mention that there are 
classically known examples of the geodesic flows on the 2-torus which 
have integrals F of degree one and two. These examples can be most 
naturally described with the use of the conformal coordinates on the 
covering plane (we refer to the book by Bolsinov, Fomenko [5] for the 
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history and discussion of this classical question) 




ds^ = A{q^,q2){dql + dql), H 
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Here A is a positive function periodic with respect to the lattice T. 
With these coordinates it can be shown that the cases of linear and 
quadratic integral for the geodesic flow correspond to those conformal 
factors A which can be written as a sum of two functions of one variable: 

A = /i(migi + ^1^2) + f2{m2qi + ^2^2) with mim2/nin2 = -1, 

where one of the functions must be constant in the case of linear in mo- 
menta integral. Such metrics are called Liouville metrics. We shall call 
a polynomial integral F„ reducible if it can be written as a polynomial 
function of the Hamiltonian H and some other polynomial integral of 
degree smaller than n, and irreducible in the other case. Let us mention 
that there are no known examples of Riemannian metrics on the 2-torus 
having irreducible integrals of degrees higher than two. This question 
is related also to the so called Birkhoff conjecture on integrable convex 
billiards in the plane ( see also [1], [I6]). In [12] (see also p3]) Kozlov 
and Denisova proved that if A is trigonometric polynomial then the geo- 
desic flow has no irreducible polynomial integrals of degree higher than 
two. Amazingly there do exist non-trivial examples of geodesic flows on 
2-sphere with integrals which are homogeneous polynomials of degrees 
3 and 4. These examples (see [S] , [S] , [S] , [B],[IS], [S]) were inspired by 
classical integrable cases of Goryachev-Chaplygin and Kovalevskaya in 
rigid body dynamics. Let us mention also that for surfaces of genus 
higher than one by Kozlov's theorem there are no nontrivial analytic 
integrals of the geodesic flow [Tl]. 




In what follows we shall work in other global coordinate system on 
the torus called semi-geodesic (or equidistant). It is built with the 
help of one regular invariant torus of the geodesic flow, call it L, which 
projects diffeomorphically to the configuration space T^. It was proved 
in [2] that such an invariant torus always exists and the Riemannian 
metric can be written in the form 



where the family {t = const} is the family of geodesies of the chosen 
invariant torus. Amazingly for our approach both coordinates will play 
a very important role while each of (c) and (s) have their own advan- 
tages. In both coordinate systems the condition of flow-invariance can 
be reduced to a quasi-linear system of equations on the coefficients of 
polynomial F„. However for the case (c) this system gets the form 





ds"^ = g^{t, x)dt^ + dx^, H 




A(f/)f/,,+A2(f/)f/,, =0 



(1) 
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while fore the case (s) it has a form of evolution equations: 

Ut + A{U)U, = ^. (2) 

Here U is a vector function of coefficients and Ai{U) and A{U) are n x n 
matrices. Let us mention here the advantage of the second system. 
Characteristics of it are always transversal to the direction while for 
the second they might rotate. This fact complicates of the analysis, 
but this is interesting open question to find suitable generalizations for 
the first system as well. 

It was proved in our recent paper ^ that the system (2) is in fact 
Rich or Semi-Hamiltonian system. Among other things this enables 
following P. Lax (see Serre's book [15]) to analyze blow up of smooth 
solutions along characteristics. Such an analysis was performed for 
other quasi-linear system (a reduction of Benney' chain) in Here 
we shall apply the same ideas to the system (2). As usual, for such type 
of systems one does not knows a-priori that the system is Hyperbolic 
or it may have regions with Complex eigenvalues. In this paper we 
shall concentrate on the case of degrees 3 and 4. Our main results are 
for these so called "elliptic" regions, we shall denote them by f2e, these 
are regions on the configuration space where the matrix A{U) has all 
eigenvalues different and two of them are complex-conjugate. Then it 
follows that for n = 3 the third one is obviously real. Also for n = 4 
the other two eigenvalues must be real as well (see later). Our main 
results for the geodesic fiow says that for n = 3, 4 the "elliptic" regions 
are "standard" that is the metric on them is classically integrable: 

Theorem 1. Let n = 3, then one has the following alternative: 

Either metric is flat in the region VL^ or is reducible on , that 
is it can he written as combination of H and Fi 

F3 = hF^ + 2k2HF^ 

for some explicit constants ki, ^2- 

Corollary 1. We have for the conformal model (c): 

Either metric p is flat on f2e or A = A(mgi + nq2) on Vt^ for some 
reals m, n; If in addition p is known to be real analytic metric on 
then A = A(mgi +nq2) everywhere on the whole torus and the flow 
p* necessarily has a first power integral on the whole torus T^. 

Theorem 2. Let n = A, then the following alternative holds: Either 
metric p is flat on Qe or F4 is reducible, that is it can be expressed on 

F4 = + 2k2HF2 + ik^H^ 

where F2 is a polynomial of degree 2 which is an integral of the geodesic 
flow on Qf, and ki are constants. 
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Corollary 2. The conformal factor A{qi,q2) can be written on in 
the form 

nil ^2 

Hqi,q2) = finiiqi + niq2) + g{m2qi + n2q2) with = -1. 

ni n2 

If in addition A is known to be real analytic then A can be written in 
such a form for all qi, q2 on T^. 

There are several main ingredients in the proof of these results. The 
first is the property of the quasi-linear system to be Rich or Semi- 
Hamiltonian. This means that it can be written in Riemann invariants 
on one hand and in the form of the conservation laws on the other 
hand. This was proved in our paper |3] . Next we use strong maximum 
principle for Riemann invariants corresponding to Complex eigenval- 
ues. Next we were able to show that for real eigenvalues the condition 
of genuine nonlinearity is satisfied. 
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2. Preparations 



In this section we explain first the geometric meaning of the charac- 
teristic polynomial of the matrix A{U) and also prove a crucial Lemma 
about it. 

The matrix of the system (2) is the following n x n (see |4j|) matrix 

A: 



A 



/ 

a„_i 
a„_i 





\ 




















Qn-l 



ai ^ 

2a2 — naQ 
3a3 — {n — l)ai 

{n - l)a„_i - 3a„_3 
nan - 2a„-2 



(3), 



where it is convenient to write 



n 



k=0 



(^k{t,x)^—rP2, 
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is the integral of degree n and the unknown vector function U = 
(ao, . . . , an-2, o-n-iY is a column vector of non-constant coefficients of 
Fni where one can show that the highest two coefficients can be nor- 
malized to be 

a„_i = g and a„ = 1. 
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The first interesting property of tliis system is tlie fact tliat its eigen- 
values have very precise geometric meaning: it can be shown that in 
order to compute eigenvalues of the system one just have to find critical 
points of Fn restricted to the circular fibre of the energy level. In other 
words let us differentiate F„ along the fibre of the energy level: 

Let Gn be the homogeneous polynomial which is the derivative of F„ 
in the direction of the fibre {H = 1/2} n T^T^. 

where the vector field v looks differently for the coordinates (c) and 
(s): In case (c) 

d d 



and in case (s) 



dpi dp2 
d Pi d 



'd(pi/g) g dp2' 

Define now usual polynomials Gn ^-nd F„ corresponding to G„ and 
Fn of the variable s in both models as follows: in case (c): 



and in case (s) 



^, Gn{s)^^Gn, Fnis) = ^Fn 
Pi Pi Pi 



Gnis) = ^Gn, Fn{s) = ^Fn. 



Pi 1 9 Pi Pi 

It turns out that with these notations eigenvalues of the matrix A{U) 
are related to the roots Sj of G„(of the model (s)) by 

\i = gsi. 

Moreover it is remarkable fact that the system (2) can be written in 
Riemann invariants: 

+ Ai(ri)^ = 0, i = 1, ..,n, 

where rj are just critical values of F„ on the circular fibres. Therefore 
for one has the following 

Notice that roots of Gn could be ±i then the Riemann invariant 
would have singularity. However this does not happen by the following 
lemma and its corollary. Notice that both statements are invariantly 
formulated but to prove it we use the model (c). 

Lemma 1. Assume Fn is divisible by H for a point on T^. Then Fn 

can he represented globally as Fn = HFn-2- Saying differently if Fn is 
irreducible integral then Fn is not divisible by H for any point. 
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Proof. This follows in fact from two identities found by Kolokoltsov in 
[TU] for conformal model (c). Write F„ and H in a complex way as 
follows ^ 

P = Pi-iP2, Fn = Y^ Aip'^-'p', H = —pp. 

It is proved in [TU] that Aq and An are holomorphic and thus must 
be constant on the torus. Therefore if F„ is divisible by H at some 
point on the torus then ^40, An must vanish at this point and hence 
everywhere on the torus Aq = An = . The claim follows. □ 

Corollary 3. For irreducible Fn, write the Gn to be derivative of Fn 
with respect to the fibre, then Gn is not divisible by H . In other words 
rki are never among the roots of Gn. 

Proof. We shall use the coordinates (c) for the proof: In complex no- 
tations p, p one has 

^ I [_dFn dFn 



2 \ dp dp 

If Gn would be divisible by pp one would have Aq = An = 0. The claim 
follows. □ 

3. Maximum principle for complex Riemann invariants 

Let us recall that we consider for the case n = 3,4 those regions on 
the configuration space where all eigenvalues are different and two of 
them are not real. In this case the third one is real for n = 3. Also for 
n = 4 the other two must be real as well. This is because any function 
on the circle must have maximum and minimum and eigenvalues as 
we explained in the previous section correspond to critical points of Fn 
on the fibre. Moreover the points of maximum and minimum must be 
different since it is impossible for Fn to be constant on a fibre. Since 
this would contradict the existence of the regular torus L (the point 
of intersection of such a fibre with L would be critical point for 
We shall call these regions "elliptic" regions. On the boundary dfle 
the complex conjugate root collide and become real. We shall use the 
strong maximum principle for the following 

Theorem 3. Let Si^2 = a ^iP be complex conjugate roots of the poly- 
nomial Gn, n = 3,4. Denote by 

Fjs) 

ri2 = u±iv = „ , ,^ 

(1 + si^yi'^ 

the corresponding Riemann invariants. Then u and v must be constants 
on VL(.. Moreover if Qe has a nontrivial boundary then f = on Qe- 

Proof. Consider first the case n = 4. Then there is no square root in 
the denominator of ri^2 and hence by the lemma of previous section u, v 
are smooth in the interior and continuous up to the boundary of Qg. 
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In the interior r = u + iv satisfies the following Vf + {ga + ig(3)rx = 
therefore 

{u + iv)t + {ga + igP){u + iv)^ = 0. 
Then denoting a — ga, $ — gP and reducing the system 

ut + aux - = 
vt + + avx = 

to 

Ut - fvt - = 

and eliminating u one arrives to the second order equation 

?)/(H/(r')/ 

Its principal part has negative discriminant thus the equation is elliptic. 
By the strong maximum principle v cannot attain maximum in the 
interior point. Therefore v must be constant and moreover to be zero 
if there is a boundary, because on the boundary v vanishes. Prom the 
equations it follows that u must be a constant as well. 

For the case n = 3, due to the square root in the formula one might 
have not a single- valued function for r. However in this case we shall 
consider instead, which is also a Riemann invariant and apply the 
same argument as above for r^. We have that r must be a constant. □ 

4. Proof of the main theorem for n 3 

In this section we prove Theorem 1 for case n=3. The proof requires 
the following: 

Theorem 4. Let Jle be region on with the property the polynomial 

Gs{s) — a2S^ + (2ai — 3)s^ + (3ao — 2a2)s — ai 

has one real and two complex conjugate roots. Then it follows that 
this region is a strip on the covering plane M.'^{t,x) with the slope A 

and on this strip the Riemannian metric admits 1-parametric group of 
symmetries g{t, x) = g{Xt — x) and therefore there exist a linear integral 
Fi=pi + \p2. 

Proof. Pirstly we have by Theorem 3 that ri, r2 are constants on Qg. On 
the domain fig ri,r2,r^ are coordinates in the space of field variables, 
ri, r2 being constant imply that third one satisfies the equation {r^)t + 
^3(^3)2: = 0, where A3 = s^g. Since ri,r2 are constant on fig then A3 
depends only on r^. We have got the simplest quasi- linear equation 

(r3)t + A3(r3)(r3)^ = 0. 
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Now the characteristic of this equations are integral curves of the vector 
field 

5^ + A3(r3(M))gj 

and r3 must be constant on these curves, hence A3 also remains constant 
along the curves. Therefore these curves are parallel straight lines. 
Take any of these straight lines which passes through interior point 
of Hg. Then it can not reach the boundary of fig ^-nd must remain 
in fie- This follows from the fact that for interior point Ai 7^ A2 and 
for the boundary Ai = A2. Along characteristic ri,i = 1,2,3 remain 
constant then also A^ because they are defined by coefficients which 
are parameterized by ri,r2,rs. So we have that fie is a strip of the 
slope A3 = const. Moreover as explained each has constant values 
along the characteristic. In particular g = g{x — X^t). This proves the 
theorem. □ 

Now we are in position to complete the proof of main theorem for 
n = 3. 

Proof. (Theorem 1). By Theorem 4 Qf, in coordinate (t,x) is a strip 
with the slope A, and At), where are functions of one 

variable. The coefficients satisfy the quasi-linear system 
(2) 

Ut + A{U)U^ = 0, [/ = (ao, ai, 02)'^, 
which for n — 3 takes the form 

/ ai 

A{U) = 02 2a2 - 3ao 
y a2 3 — 2ai 

Since U — U{x — Xt) is in the form of the simple wave then U' is 
A-eigenvector oi A{U). Then one has 

~ AOg 

a2aQ + 20202 — 3aoa2 = Xa'^ (4) 
a2fl'i + (3 — 2ai)a2 = Aa'g. 

These differential equations can be solved as follows: 
Divide the last equation of (4) by ol have 

ai 3 — A 1 

^2 2 02 

and so 

ai = (3-A)/2 + cia2. (5) 
Divide the second equation of the system (4) by 03 to have 

(ao/al)' = (1/02)' + '^«i/«2 
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This expression together with (5) yields 

ao 1 — Aci 

which means that 



+ C2, 



ao = a2(l - Aci) + C2a2. (6) 

On the other hand substituting a\ from (5) into the first equation of 
the system (4) one gets 

, Ci o 3 — A , . 

Aao = —a^ H —a^ + Cg. (7) 

Ehminating ao from the equations (6), (7) one gets certain third power 
polynomial on a2 which vanishes. Then there are two possibilities: ei- 
ther the function a2 is a constant and then the metric is flat (remember 
0^1 = 9) or coefficients of this polynomial must vanish. But this yields 
the identities: 

-n _ 3(A-1) _A-1 
Using them one can easily verify the following explicit identity 

I. I 3-A 

^^ = '^ = ^' ^■^ = ^ 
(where F3 = + ^^pjp2 + PipI +pI)- 

Let us remark here that the case A = means a, = ai{x) in particular 
g — g{x) but then p — g'^{x)dt^ + dx^ is obviously flat metric. □ 

Proof. (Corollary 1). Let us remark first that the fact that the cubic 
integral can be explicitly expressed through the first power integral is 
absolutely necessary for the proof. This is because the elliptic domain 
could be proper subset of the torus. In such a case it is not clear why 
its coefficients should be constants. Therefore we proceed indirectly as 
follows. 

By the previous theorem we have that the metric p posses linear in 
momenta integral Fi on the set f2e and moreover F3 can be expressed 
through Fi and H. Write Fi = bo{x, y)pi+bi{x, y)p2- Using the identity 
F3 = kiFf + 2k2HFi and Kolokoltsov constants for F3 we get bo and 
bi must be constants (obviously at least one of them is not zero). But 
then the last equation of quasi-linear system of the conformal model 
gives: 



therefore A = A{biqi — 6092)- If we know that A is of this form on an 
open subset of and if A is analytic then obviously A is of this form 
on the whole torus T^. This completes the proof. 
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□ 



5. The case 



We shall split our proof in the two cases: 

The first is the case when either has a nontrivial boundary or the 

Riemann invariants ri,2 becomes real at some point on the torus. And 
the second case (in the next section), when coincides with the whole 
torus and ri 2 are not reals everywhere on it. 

Theorem 5. Let F4 : T*T^ M. be a polynomial of degree 4 such that 
{F4, H} = 0. Denote by fig the domain on where polynomial G4 has 
two complex-conjugate and two real distinct roots. Then if ft^ 7^ or 
Jig = but one knows that Imri^2 = then F4 can be expressed on Qg 
as follows 



where F2 is a polynomial of degree 2 which is an integral of the geodesic 
flow. 

Proof. Denote by Si^2 = « ± «/9 the pair of complex-conjugate roots 
of G4, and by ri^2 = u ±iv the corresponding Riemann invariants. It 
follows from the Theorem 3 that v = and m is a constant on Q^. So 
we have that ri^2 is a real constant. Denote it by r. Then we claim 
that F4 - 4ri7^ can be factorized F4 - 4rif^ = KM where K and M 
are real polynomials of degree 2 and 



must be divisible by K. Therefore there are two possibilities: 

1. M and K are relatively prime. In this case LyK must be divis- 
ible by K. Notice that for any quadratic polynomial K, L^K always 
have real roots. This is because any function on the circle must have 
minimum and maximum. Thus the only possibility in the first case is 
LyK = 0. But this means K is proportional to H. One may assume 
K = H (by taking the coefficient of proportionality to M). So we have 
got F4 — ArH'^ — HM, therefore M must be integral of the geodesic 
flow degree 2. 

2. In this case M is divisible by K, i.e. proportional to K, M — cK 
(where c{t,x) is a function). So we have F4 — ArH'^ — cK'^. One has 
equating the coefficients the coefficients of (^3) at both sides: 1—r — c. 
So c is a constant and we have 



F4 = + 2k2HF2 + AksH^ 




Since {a ± i(3) are roots of G4 then 



G. 



'4 = LyF = {LyK)M + K{LyM) 



F4 = ArH^ + (1 - r)K^. 
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Again K must be an integral of degree 2. So in both cases we proved 
that F4 is reducible. This proves Theorem 5. □ 

By the same method we can prove the following 

Corollary 4. Under the conditions of Theorem 5 the conformal factor 
A(gi,g2) can be written on fie in the form 

HQi^Q2) = fii'miqi+niq2)+f2im2qi+n2q2) with mima/ning = -1. 

Moreover, if A is known to be real analytic then A can be written in 
such a form for all (gi, ^2) on T^. 

Proof. We bypass the difficulty exactly as in the previous corollary. As 
we proved one knows the existence of integral F2 on the domain Q^, 
such that F4 can be expressed as function of H and F2, then write F2 
in conformal model (c) in the complex form F2 = h^p^ + hipp + &2P^- It 
follows from Theorem 5 and Kolokoltsov identities that 60? ^2 must be 
constants on fie, = A + iB, b2 = A — iB and bi is a real function. 
Then the equations on bi and A are the following: 

{biA)g, = -2AAg, - 2BAg„ (61A),, = -2BAg, + 2AAg,. 

Eliminating (&iA) we get 

B^qiQl ~ 2^Aq^g2 — BAq^q,^ = 

which imply the result. □ 



6. The case = and ri^2 are not real 

This is the most difficult case which is left. In this case we shall 
prove that each one of the real eigenvalues is genuinely nonlinear in 
the sense of Lax. This fact, together with the property of our system 
to be Rich (or Semi-Hamiltonian), will enable us to establish blow-up 
of the derivative Vx unless it vanishes. Such a proof was first given in 
[3] for other system. 

Theorem 6. Assume Qe = T^, d'^d assume that for all {t,x) the poly- 
nomial G4 has 4 distinct roots, 2-complex conjugate Si^2 = OL±i(i and 
2 real S3.4. Assume also that the imaginary part of Riemann invariants 
ri,2 does not vanish. Then the real eigenvalues X^^^ = (753,4 are nec- 
essarily genuinely non-linear and therefore the corresponding Riemann 
invariants are constants. In particular all ai must be constant, and so 
the metric is flat. 

We proceed as follows. Let us first subtract from F4 — AH^ in order 
to make the coefficient of P2 vanish, all the other we shall denote by 
the same letters. 
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Lemma 2. Let A3 = gs^ be an eigenvalue of our quasi-linear system 
then the derivative |^ can be computed 

^ ^ ~ o/^,r , (("3 + oijsg + 4aoS3 - (oi + 03)) 
era 5^^4(53) 

5" = (S3 - Si)(s3 - S2)(S3 - S4), 

^4(5) = 035^ + 2a2S^ + 3(ai — 03)5^ + (4ao — 202)5 — oi. 
Proof. We have 

SAs 19(5303) 9a3 / aS3 Sfli 



drs dr^ dr^ \ ^ dai dr^ 



Since S3 is a root of G4 = then 

ds^ _ 4s3 as3 _ 1 - 3si 
^ ~ ~5S) ' ^ ~ 6^(53) ' 
9^3 ^ 2s3(l - gj) as3 _ sl{sl - 3) 
5a2 G'^{s^) ' 9a3 ^^(53) 
In order to find ^ we need to calculate |^ : 

ar-3 dr-3 



9r j / \ dttj 
Prom the previous identity we obtain 

dttQ _ (1 + s|)^SiS2S4 9ai _ (1 + s|)^(siS2 + + S2S4) 

da2 _ (l + gi)^(gl+g2 + g4) _ (l_fs|)^ 

9r3 S" ' dr^ S 

The fact that s^, i = 1, . . . , 4 are roots of G4 gives us 

2a2 

Sl + S2 + S4 = S3, 

as 

3(01 + 03) , . 

S1S2 + S1S4 + S2S4 = 53(01 + S2 + Si), 

03 

3(«i + «3) /2a2 , \ 

H S3 K S3 , 

03 V «3 / 



S1S2S4 = 

0353 

and 



dttQ _ (1 + slfai dai _ (1 + s^)^(3(ai - 03) + 202^3 + Qssj) 
drs a^SsS ' 9r3 036" 

da2 _ (l + ,sg)^(2a2 + »3.S3) 

9r3 035* 
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From here we have 

^ = -—^ — ^(2a34+4a2S3+3(ai-3a3)s^-4(ao+a2)s3+ai+3a3). 

Cra 5*^4(53) 

The identity 

2G4(s3) - (20383 + 402^3 + - 303)53 - 4(ao + 02)53 + Oi + 303) 

(03 + 01)53 + 4oo53 - (oi + O3). 
proves the Lemma 2. □ 

Genuine nonhnearity is proved in the following lemma which is in 
fact general fact about critical values of the polynomials on the circle. 

Lemma 3. Let F4 be a polynomial of degree 4 Fi = Xli=i ( ^ ) P'^ 

considered on the circle H — ^ (^(^) ~ ^' ^^''^ote by G4 the 

derivative of F4. Let Si be the roots of such that 51,2 are complex 
conjugate and 53^4 are real distinct. Let ri denote critical values, with 
ri,2 are complex conjugate (not real) and r^^^ real. Then it follows that 
for Xi = gsi 

/ and / 0. 

ors or^ 

Proof. In order to prove the lemma denote 

51,2 = a ±i^, P¥^0. 

We have to show that in such a case polynomials (74(5) and 

7 = (03 + ai)s^ + 4oo5 - (oi + 03) 

are relatively prime (notice that 7 has only real roots). This can be 
done as follows. Divide G4 by 7 with a reminder R. If G4 and 7 have 
common root then R has to be of degree zero or to be equal to zero. 
Dividing explicitly one gets for it! the following expression: 

R 

2(oi + O1O3 - oi (40002 + o|) + 03(800 - 4oo02 - ai))(ai + 03 - 4oo5) 

(oi + 03)3 

Notice that R can be a number in two cases only: 
Case 1. R ^ but oq = 0. 

In this case 7 has two roots ±1 and the value of G4 is the same for 
both of them 

G4(±l) = 2(01-03). 

But this means that if one of them is a common root of G4 and 7 then in 
fact both of them are. But in the first case R is not zero, contradiction. 
Case 2. = 0. 
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In this case G4 is divisible by 7. Denote a ± —l/n the roots 
of G4 where the last two are the roots of 7. With these notations one 
obtain from the Viete's formula the relation between them 

+ — 1 = a(fi — 1/ fi). 

Notice that if a = then /3 = 1 and this case is excluded by Lemma 1. 
Moreover the Viete's formulas together with this relation lead to the 
following expressions 

as Aa ' 

— — a + p , 
03 

02 1 - - /52 



03 2a 

Using these formulas one can substitute them into the value of Rie- 
mann invariant of the root a + if3 to get by direct calculations that its 
imaginary part vanishes identically: 



Im- 



F4{a + i(3) 
l + (a + «/3)2)2 



oo + 01(0; + + 0,2(0 + + 03(0; + _ n. 

But this is not possible by the assumptions. 

The exceptional cases can be treated analogously: 
Case 3. It could happen that 7 is identically zero, i.e. 

Oi + 03 = 0, Oq = 0. 

In this case — = s'^ + 2ks^ — &s^ — 2ks — l, where k = —. One can check 

as ' as 

that this is impossible because such a polynomial has 4 real roots. One 
can see this for instance through Q^^g2^'2 which in this case is 

F4 + (1 2^' + 



(1 + S2)2 (l + s2)2 

But the expression in the nominator has zero as a root and two more 
real roots. Therefore G4 must have 4 real roots. Contradiction. 

Case 4. In this case it could happen that 01 + 03 = and the common 
root of G4 and 7 is s = 0. But then ai = 0. Then also 03 = but this 
is not possible since = g is always positive. This finishes the proof 
of lemma. 

□ 



CUBIC AND QUARTIC INTEGRALS 
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Proof. (Theorem 6). It follows from Lemma 2 that real eigenvalues are 
genuinely non-linear. It was shown in our previous paper [1] that our 
system is rich (Semi-Hamiltonian). This property is crucial because it 
enables to use Lax analysis of blow up along characteristics as it was 
done for Benney chain in [3j. This method uses Ricatti type equation 
which (ra)^; and (r4)x must satisfy, and if one knows that A3, A4 are 
genuinely non-linear (Lemma 2) then {rs)^ = [r^jx = since otherwise 
there is a blow up after a finite time for Ricatti equation. We had 
before already ri = const, = const. Thus all Riemann invariants are 
constant and so aj. In particular g is constant and the Riemann metric 
is flat. □ 



7. Concluding remarks and questions 

1. It would be very interesting to know if our Semi-Hamiltonian 
system (2) is in fact Hamiltonian and to find Dubrovin-Novikov bracket 
of hydrodynamic type (see [7]). 

2. In this paper we show that in the case n = 3, 4 the system (2) 
standard in the elliptic domain Q^- It follows from our main theorems 
that in the analytic case we can assume that for n = 3, 4 hyperbolic do- 
main is whole torus T^. It follows from [T7] that in Hyperbolic domain 
Semi-Hamiltonian systems can be integrated by Tsarev's generalized 
hodograph method. It would be very interesting to apply this method 
to the system (2). 

3. It would be natural to find the blow up mechanism for the quasi- 
linear system in the conformal model (c). Technically this is not obvi- 
ous because the characteristic fields may rotate on the torus. 
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